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Abstract 

It was predicted that local spin polarization in a ring of five dipolar coupled 
spins should present a particular fingerprint of quantum interferences reflect- 
ing both the discrete and finite nature of the system [Phys. Rev. Lett. 75 
(1995) 4310 ]. We report its observation for the proton system of a (CsHs^Fe 
molecule using a rare 13 C as local probe. Novel high frequency (~ 60/cHz) 
polarization oscillations appear because incomplete 13 C- 1 H cross-polarization 
transfer splits the polarization state, in a portion that wanders in the proton 
system and one that remains in the 13 C. They interfere with each other after 
rejoining. 
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I. INTRODUCTION. 



Mesoscopic systems |I| constitute a very active topic because their transport properties 
manifest quantum mechanics in its full splendor and they open new possibilities for the design 
of devices and finely tailored systems which started with the tunnel diode . These involve 
interference processes both in the space and time domains. The last case includes the conse- 
quences of a finite tunneling time in resonant devices P and the propagation of a localized 
excitation in a finite random media. Here, a simple average over impurities configurations 
produces an apparent diffusion, however, a non-linear cx-model calculation [Q predicts strong 
long time correlations, meaning that a local excitation returns as a mesoscopic dynamical 
echo after diffusing away and reaching the system boundaries. The observation of these 
effects in the time domain is difficult because excitations propagate through the system in 
a time scale that is in the frontier of present technological possibilities. However, a similar 
phenomenon was also predicted || for a local excitation in a system of nuclear spins with 
magnetic dipolar interaction, which evolves with characteristic times on the hundreds of 
micro-seconds and can be tested by the powerful NMR spectroscopy. While the discrete 
nature of spin systems leads to a fingerprint of quantum interferences in the spin dynamics 
of both finite and infinite systems, the survival of mesoscopic phenomena might not seem 
obvious since, at high temperature, different configurations must be averaged and it is often 
described as a diffusive process. Therefore, in order to illustrate the physics we resort to an 
XY model, where exact analytical results can be obtained (technical details sketched bellow). 
Fig. 1 shows the evolution of the local magnetic polarization M 7 (t) for a ring of 21 spins and 
an infinite chain at high temperature. At short times a series of quantum beats, associated 
with the discrete nature of the system, are developed. At around 2300/is it appears a revival 
of the polarization. A study of polarization at other sites shows that the wave packet splits 
in two parts winding around the ring in opposite directions and the maximum appears when 
they meet again in the original site. Similarly, reflections at the edges of linear chains pro- 
duce well defined maxima. While we called this phenomenon quantum dynamical echo 
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in analogy with the prediction for quantum dots ||], to avoid confusions that might arise 
when applied to NMR experiments, we will use mesoscopic beat henceforth. In contrast 
with ordinary quantum beats, its mesoscopic origin is shown by a time scale proportional 
to the system size (see inset of Fig. 1). In dipolar and J-coupled systems the range of sizes 
that present mesoscopic beats is much more restricted than what the XY example suggests. 
For very small systems they become indistinguishable from quantum beats, while for large 
systems they are attenuated progressively until they merge in the background of the other 
beats. Our finding that local polarization in a system of five dipolar coupled spins arranged 
in a ring configuration should retain quantum beats, was surprising since experiments |7j || 
in such a system did not show any signature of these quantum phenomena. This motivated 
our present experiments under precisely controlled conditions. 



The NMR radio-frequency pulse sequence developed by Zhang, Meier and Ernst |7j 
(ZME), involves cross-polarization (CP) ||, to achieve a two-way transfer of polarization 
|T0|j among the spin S of a rare 13 C (1.1% abundance) and the I\ spin of a directly bonded 
1 H . Then, the 13 C can be used as a local probe that injects magnetization in the proton and 
later captures it. The set of coupled protons within a molecule constitutes the mesoscopic 
system where the spin dynamics can be monitored, while the rest of the crystal constitutes 
a 'weakly' interacting reservoir. We measure this spin dynamics in a polycrystalline sample 
of ferrocene, (CsHs^Fe, (Fig.2). Since the dipolar interaction depends on the angle between 
the internuclear vector and the magnetic field, in a general molecule each pair of nuclei can 
have a different interaction parameter. Therefore, the resulting Hamiltonian is a particular 
case of: 



j>k k 

where subscripts indicate spins. For a = 0, it defines an XY model, a = — \ describes the J- 
coupling or Heisenberg model, and a = 1 is the truncated dipolar Hamiltonian. In ferrocene 



II. THEORY AND EXPERIMENTAL METHODS. 




(1) 
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at room temperature, the rings perform very fast rotations around the C5 molecular axis 
(tr w 1(T 12 s) |TI|. In this last case djk are the time-averaged ]I2| interactions parameters 
depending only upon the angle between the rotation axis and the magnetic field. For each 
initial state, | z),with site 1 polarized, the probability of finding the same site polarized in 
the state (/ | after a time t is: 

; 2 



(f I exp[--H H t] I i) 



(2) 



from which the polarization can be calculated summing over all the Ni and Nf possible 
initial and final states: 

~N f Ni 



M\t) 
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This is the magnitude which presents quantum mesoscopic beats as function of t. For the case 



a=0, ( XY-model) an exact mapping Jl3 to a non-interacting fermions system, allows to sum 

2 

up the terms in Eq. (|3|) and write the magnetization just as M (t) = (1 | exp[— VHnt) | 1) , 
the square modulus of a single particle wave function, where state | 1) has spin I\ "up" and 
all the others " down" . Even when the a^O cases map to interacting particles, the essence 
of this excitation dynamics is retained until intermediate times. 

The best condition to observe quantum beats || occurs when the polarization is quan- 
tized along the external magnetic field (laboratory frame). This is because the neglect of 
non-secular terms in the dipolar Hamiltonian is better justified and, since dipolar interac- 
tions are maximized, the time scale of quantum dynamics shrinks minimizing the effects 
of other interactions leading to relaxation. The complete pulse sequence is schematized 
in Fig. 3. A (ir/2) x pulse on the abundant X H spins system creates a polarization that 
is transferred during tc to the rare 13 C system when both are irradiated at their respec- 
tive resonant frequencies with field strengths fulfilling the Hartmann-Hahn |9j condition 
jiBfj = 7s-Bis = ^is (here u> is = 2tc x 44.6kHz). After decay of the proton spin coherence 
during time t$ = 1ms the most relevant part begins: A) The magnetization from an initially 
y polarized 5* spin is transferred to the y-axis of the 7 X spin through a CP pulse of duration 



td = 85/is (the shortest maximizing the polarization transfer |TJj for the selected orienta- 
tion). B) A (tt/2) pulse tilts the polarization to the laboratory frame. C) The /-spins 
evolve freely during a time t 2 - Thus, the relevant evolution of the spin-diffusion sequence 
occurs in the laboratory frame while S— irradiation prevents the system-probe coupling. D) 
A (vr/2)_ a , pulse tilts the polarization back to the rotating xy plane. E) Another CP pulse 
of length t p = td is applied to transfer back the polarization to the x-axis of S. F) The S 
polarization is detected while the /-system is kept irradiated (high-resolution condition) . An 
important feature in the pulse sequence of fig 2 is that while during t 2 the spins evolve with 
Hamiltonian ([J), during CP periods td and t p the proton spins evolve with a Hamiltonian 
which after truncations is TL' n = — [^\TC U . Thus a part of the free evolution time t 2 (around 
— [-^(td + tp) =42.5/xs) is spent going backwards in time (the Loschmidt's daemon [TJj 
is working!) compensating the undesired evolution. In the ZME experiment an additional 
polarization evolution in the rotating frame is allowed by keeping the r.f. in the proton 
system during a time ti(i.e. after CP shown in A and before the tilting pulse B). Therefore 
a time t 2 = |ti + t m is required to reach the local polarization maximum (polarization echo), 
after which it evolves according to Eq. 3. 



III. RESULTS AND DISCUSSION. 

In a polycrystal all the angles between the molecular rotation axis and the magnetic field 
are equally present. Therefore the 13 C spectrum has a shape typical of axially symmetric 
systems with a very well developed peak corresponding to molecules with their rotational 
axis lying in the plane perpendicular to the field. Then, if we select the signal intensity at this 
exact frequency, we are monitoring the spin dynamics only on those molecules which have 
the same intramolecular interactions (dip = 1576Hz x 2irh). Besides, the carrier frequency 
is set to have this peak exactly on-resonance avoiding time evolution with chemical shift. 
The experimental results with the sequence of Fig 3 are shown in Figure 4 with circles. 
Diamonds correspond to a run with the full sequence [7J of ZME with t% = 80/is. It is 
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presented here as a function of ti with a shift — ti/2 and properly normalized. The line 
is obtained by a spline fitting. For short times we see that superimposed to a parabolic 
curve with a maximum (polarization echo) around t m there are high frequency oscillations. 
This is a novel quantum phenomenon that we discuss bellow. Local polarization decreases 
until a quantum beat appears with a clearly developed maximum at ~370/is. A second 
maximum which, on the basis of numerical analysis, we identify with the mesoscopic beat, 
develops at 520/is but it is overcome by an overall attenuation with a characteristic time 
of about 500/is. This attenuation contains both the interaction with neighbor molecules 
and decoherent processes. The inset shows a sequence of ideal (without couplings to 13 C) 
evolutions, M (t), calculated with Hamiltonian (|]) for a single ring and complete molecule 
(10 spins in Fig. 2) as described in ||. The single ring is the upper dotted curve and is the 
case where our previous study of different correlation functions allowed the identification of 
the second peak at 580/zs as a mesoscopic beat ||. For the molecule we distinguish four 
situations: rings rotating independently but keeping a) a staggered configuration , or b) an 
eclipsed one ; molecule rotating rigidly with their rings in c) a staggered configuration or d) 
an eclipsed one. Each of these gives different sets of time averaged interaction parameters and 
represents a progressive increase in the inter-ring interactions which produces the consequent 
attenuation of the interference phenomena. An important conclusion to be drawn from this 
sequence is that while changes in the short time evolution, M 7 (t) ~ 1 — |Aco> 2 x t 2 , are not 
much apparent, the effect of those interactions in blurring out the long range interferences is 
decisive. This is evident in the decrease of the mesoscopic beat. The above arguments also 
lead to interpret the additional attenuation shown up in the experimental curve as due to 
intermolecular interactions. Therefore, although the experimental data can not be fitted to 
our simulation with small number of spins, the comparison of these theoretical curves with 
the experimental one establishes that the rings rotate independently, a result compatible 
with conclusions drawn from proton T\ measurements ||16|| , but can not distinguish between 
staggered and eclipsed configurations. The clear correlation between the positions of the 
experimental and calculated peaks indicates that this is the first observation of a mesoscopic 
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beat in a system of spins with magnetic dipolar interactions. 

The short time oscillations can be explained as follows: Polarization initially at the 13 C 
nuclei has to be transferred to the proton system. Since systems with more than one proton 
have an energy uncertainty due to the inter-proton interaction, it is not possible to obtain an 
'exact' Hartmann-Hahn condition flT0 |. Therefore, the initial 5* polarized state splits in two: 



Polarization stays at S ( state A) with a probability amplitude a = cos <fi, or it jumps to the 
proton system (state B) with a probability amplitude b = i sin 0. During the free evolution 
time the state A acquires an additional phase factor due to the spin lock field. When a new 
CP pulse allows these states to rejoin, the phase produces the oscillation in the measured 
intensity. The simplest model showing this phenomenon considers a situation in which the 
incomplete (\b\ 2 < 1) CP transfer occurs instantaneously, with both polarization transfers 
along the y-direction, and the spin dynamics of only two 1 H with an XY interaction. Then 
the polarization injected in one of them at time ti = will be find again there with a 
probability amplitude ms(t) = cos[|f2 t], with Q obtained as the difference between singlet 
and triplet eigenenergies. The ideal magnetization is then M 7 (t) = |ms(t)| 2 . Because 
of r.f. irradiation during t%, the A polarization amplitude is = exp [iuis t] . The 

magnetization measured at the 13 C non-ideal probe after a new contact pulse can be written, 
after some algebra, as 

M s (t) = |6| 4 M 7 (t) + 2 x 6 2 cos[ifi t] x a 2 cos[a; ls t] + |a| 4 . (4) 

Here we identify an interference term proportional to the waiting amplitude, mg(i). To 
emphasize the effect we re-write this as: M s {t) = |6 2 m^(t) + a 2 m^(t)| 2 . The first term 
represents the probability amplitude for the S polarization to be transferred to the proton 
system, times the probability amplitude for the magnetization to remain in the original 
proton site (waiting amplitude), times the probability amplitude to be transferred back to 
S; the other is the probability amplitude for the polarization to remain in S, evolve and still 
remain in the 13 C. In the actual experiment (Fig. 3) the second transfer is to the x-axis, 
this cancels out the last term in Eq. @) and decreases by a factor of 2 the interference 



term. That gives a coefficient \a\ 2 / \b\ 2 for the relative importance of the oscillations. From 
these considerations, the experimental oscillation amplitude implies that there is about 90% 
efficiency in the polarization transfer from 13 C. This is also consistent with the experimental 
lower bound of 85% efficiency for each CP transfer obtained from the ratio of intensities 
measured with ti — and with a sequence in which portions A-E have been suppressed. 
This simple case has the virtue to present a structure that is maintained for more complex 
proton systems: the first and third terms in Eq. (gj) represent the classical effect of sum of 
probabilities, the second term is an interference one. In general, the frequencies involved 
in the first two terms are not the same. Those contributing to the first have the form 
Qij = (Ef — E^)/ti, where the superscript in the eigenenergies of Hamiltonian (H) indicates 
the number of "up" spins in the proton system, while the interference term contains the 
frequencies Vt'^ = {E^ +1 —E^)/fi i.e. the subspaces mixed by CP. Hence the expression of the 
observed polarization M s (t2) is not simple. A better model takes into account the evolution 
in the proton system during the CP periods. Its numerical integration gives, in agreement 
with the experiments, both an overall shift in the oscillation toward frequencies higher than 
(jj\s and an oscillation maximum which precedes that of M 1 ^ — t m ). Further comparisons 
between models and experiments in simple systems should provide useful information about 
decoherence processes and efficiency of cross polarization pulses. On a more speculative 
scope, we notice that the described manipulation of polarization is of the type proposed 
to implement simple quantum computations |l7j and its further understanding could also 
contribute to this developing field. 

One might wonder why these wealth of quantum phenomena were not seen in the ZME 
experiments performed in a single crystal 0. In order to observe the interference between 
polarization pathways it is important to take data at short enough time intervals. The 
experiments of ZME were performed at about the characteristic period (2ir/uis) and that 
explains why they could not see this phenomenon. Note that use of higher r.f. power would 
make the observation even more elusive. Finally, as mentioned above, in a polycrystal we 
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are observing simultaneously all the molecules which have their rotational axis in the plane 
perpendicular to the magnetic field; among them there are some molecules that have the 
interactions with their neighbors minimized while other that have them maximized (in a 
single crystal this shows up as an orientation dependent second moment Au 2 ). While rings 
with stronger interaction contribute with polarization evolutions which do not have quantum 
beats but decay steadily much as the plots c) and d) in the inset of Fig. (3), the molecules 
in which the characteristic inter-ring vectors form an angle close to the magic angle can 
contribute with well defined quantum beats. The importance of the interaction of neighbor 
molecules is appreciated in the overall decay of polarization observed in the experimental 
data in contrast with the finite asymptotic polarization in the exact solution of an isolated 
molecule. While our experimental data are a superposition of these typical behaviors, the 
reported single crystal data, are probably arising from an orientation favoring intermolecular 
interactions. Therefore we expect that by choosing the appropriate orientation in the single 
crystal a mesoscopic beat better developed than the one shown here should be observed. 
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Figure 1: Evolution of the local magnetic polarization M 7 (t) for a ring of N =21 spins 
at infinite temperature with an XY (a = 0) interaction d = 1576 Hz x 2nh. Fine line is the 
square modulus of the J Bessel function, exact solution of an infinite chain. At short times 
a series of quantum beats (at around t n ~ z n h/d , z n the zeroes of Ji) are developed which 
decay with a t^ 1 law, indicating the excitation propagation. At around 2330/is = t m .b.. it 
appears a revival of the polarization corresponding to the excitation winding around the 
ring. The inset shows the mesoscopic beat time, t m .b., for different ring sizes. Notice that 
t m .b. oc N x a/v, where a is a lattice constant and v ~ ad/% an effective group velocity. 

Figure 2: Ferrocene molecule in a staggered configuration. Fe atom is at center. Hy- 
drogen atoms are labeled starting from the one which has a rare 13 C directly bonded. The 
eclipsed configuration is obtained rotating any ring 36° around the C5 axis. 

Figure 3: Pulse sequence for proton spin 'diffusion' in the laboratory frame. A (tt/2) x 
pulse on the abundant / spins creates polarization. This is transferred during tc to a rare 
13 C spin S. After decay of proton spin coherence during ts- A) The S magnetization is 
transferred to the y-axis of I x through a CP pulse of duration t d . B) A (n/2) x pulse tilts 
the polarization to the laboratory frame where C) it evolves freely during t 2 - D)A (ir/2)_ x 
pulse tilts it back to y and E) a CP pulse of length t p = td transfers the polarization to 
the x-axis of S where F) it is recorded. Dashed and dotted lines show main and secondary 
probability amplitude pathways for the polarization. 

Figure 4: Evolution of a laboratory frame X H spin polarization in ferrocene as detected 
in a directly bonded 13 C spin. Notice the maxima at ft;370/zs and ~ 520/is (mesoscopic beat) 
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and quantum interference oscillations of period 16/is. Inset: Calculated ideal evolution 
of local 1 H polarization in an isolated ring (dots) and in the complete molecule (see Fig 
2) where the rings rotate independently in: a) staggered (full line) or b) eclipsed (small 
dashed) configurations, and the rigid rotating molecule keeping: c) staggered (long-dashed) 
or d) eclipsed (dot-dashed) rings. 
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